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Abstract. We present new calculations for transition probabilities of Fe IV, with much more extensive datasets than heretofore 
available. The large-scale close coupling R-matrix calculations yield 1,798 LS bound states with n < 11 and / < 9, and corre- 
sponding transition probabilities for 138,121 dipole allowed transitions in the form of oscillator strengths /, line strengths 5, 
and A-coeffficients for a variety of applications. This represents the largest R-matrix dataset in LS coupling for any ion under 
either the Opacity Project or the Iron Project. Through algebraic transformation of the LS multiplets, a total of 712,120 dipole 
allowed fine structure transitions for Fe IV are obtained. Observed transition energies are used together with the energy indepen- 
dent line strengths to derive the /- and the A-values; the adopted algorithm used calculated energies for the remainder. Present 
results show significantly better accuracy for the important low-lying states than previous calculations. Monochromatic and 
mean opacities for Fe IV are computed and compared with those obtained using the Opacity Project data. We find differences 
which could have important consequences for several astrophysical applications involving low ionization stages of iron. 
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1. Introduction 

It is well known that low ionization stage of iron are of great 
interest for astrophysical spectroscopy and models of stellar at- 
mospheres (Hubeny etal. 2003), novae (McKenna etal. 1997), 
photoionized H II regions (Dopita and Sutherland 2003, Rubin 
et al. 1997), and active galactic nuclei (Sigut and Pradhan 2003, 
Sigut etal. 2004). The ion Fe IV is of particular interest since 
its ionization potential is ~ 54.8 eV, nearly equal that of He II 
at 54.4 eV. Thus, for example, in a photoionized H II region 
such as the great diffuse nebula in Orion, stellar ionizing flux is 
effectively shielded by Helium to prevent extensive ionization 
of Fe ions beyond Fe IV. Whereas hot central stars of high- 
excitation planetary nebulae may well ionize beyond Fe IV 
(e.g. Hyung and Aller 1998), ordinary main sequence O, B stars 
have a 54 eV 'cut-off' that attenuates the ionizing flux. Another 
example is the anomalous Fe II emission from active galactic 
nuclei. Non-local-Thermodynamic -Equilibrium (NLTE) spec- 
tral models need to consider the first four ionization stages, 
Fe I - IV, with Fe IV as the dominate stage in the fully ionized 
zone (e.g. Sigut etal. 2004 and references therein). In previous 
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studies at Ohio State under the Iron Project, we have reported 
a large amount of radiative and collisional data for the ioniza- 
tion stages Fe I- VI, including fine structure transition probabil- 
ities for Fe II and Fe III (references at: www.astronomy.ohio- 
state.edu/~pradhan). The aim of this paper is to report similarly 
new calculations with improved accuracy for Fe IV including 
fine structure that is often required for spectral diagnostics and 
modeling. 

However, study of Fe IV is difficult and complexities arise 
due to its electronic structure with 5 electrons in the half-filled 
open 3d shell (ground configuration ls 2 2s 2 2/? 6 3.s 2 3/7 6 3af 5 ); 
electron correlation effects for the many excited states are diffi- 
cult to represent. Prior to the Opacity Project (OP, 1995, 1996), 
investigations on radiative processes for Fe IV were carried out 
through atomic structure calculations. These include radiative 
decay rates for forbidden transitions by Garstang (1958), which 
are available from the evaluated compilation by the National 
Institute for Standards and Technology (NIST: www.nist.gov) 
and Rubin and Froese-Fisher (1998). Fawcett (1989) has com- 
puted the data for many allowed transitions using a semi- 
empirical method due to Cowan. 

Ab initio calculatons in the close coupling approxima- 
tion using R-matrix method for the radiative processes were 
carried out under the Opacity Project (OP, 1995, 1996) by 
Sawey and Berrington (1992), and later under the Iron Project 
(IP, Hummer et al. 1993) by Bautista and Pradhan (1997). 
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Further R-matrix calculations were reported for the total and 
state-specific electron-ion recombination rate coefficients using 
the unified treatment by Nahar, Bautista, and Pradhan (1998). 
Sawey and Berrington (1992) produced radiative data for pho- 
toionization and oscillator strengths for three symmetries: sex- 
tets, quartets and doublets; these data are accessible through 
the OP database, TOPbase (Cunto etal. 1993). Later work of 
Bautista and Pradhan (1997) took more accurate account of 
electron-electron correlation effects, showing important reso- 
nant features in cross sections missing in the former investiga- 
tion by Sawey and Berrington (1992). However, owing to com- 
putational constraints, Bautista and Pradhan (1997) entirely 
omitted the doublet symmetry which has the maximum num- 
ber of bound states. All previous R-matrix results are in LS 
coupling and consider only the dipole allowed LS multiplets; 
no fine structure splitting has been taken into account. 

Present work reports a more complete set of oscillator 
strengths and radiative decay rates of Fe IV for both the 
LS multiplets and allowed fine structure transition. Similar 
results have been presented earlier for Fe II (Nahar 1994), 
Fe III (Nahar and Pradhan 1996), and for some other ions. 
Present calculations consider correlation effects accurately, as 
in Bautista and Pradhan (1997), and bound states of all sym- 
metries, sextets, quartets, and doublets. At the present time it 
is computationally prohibitive to carry out full Breit-Pauli R- 
matrix calculations for fine structure transition probabilities in- 
cluding relativistic effects in an ab intio manner. Among the 
low ionization stages or iron, Fe I - V, there has been only one 
such previous effort, for Fe V (Nahar and Pradhan 2000, Nahar 
et al. 2000). challenging computationally. While efforts are un- 
der way to resolve the computationally challenging problems 
in relativistic calculations, which are about an order of magni- 
tude more intensive, there is considerable need for fine struc- 
ture data. We employ an algebraic recoupling algorithm to ob- 
tain the fine structure transition probabilities for Fe IV, using 
observed energies where available, as in previous IP works on 
Fe II and Fe III. 

Finally, we employ the new radiative dataset to compute 
monochromatic and mean opacities for Fe IV, both as a test and 
comparison with the earlier OP work, as well as to investigate 
possible effect on the astrophysically important problem of iron 
opacities. 

2. Theory 

The R-matrix method as applied to the OP and IP work enables 
the consideration of many excited bound states and transitions, 
typically for all states with n < 10. General details of the pro- 
cedures employed have been described in earlier papers, such 
as for Fe II (Nahar 1993), for obtaining oscillator strengths (/) 
and radiative decay rates (A) for dipole allowed fine structure 
transitions (AJ = 0,±1 with AS = 0, even parity <-> odd parity) 
by algebraic tansformation of LS multiplets. The methodology 
of the close coupling (CC) approximation employing the R- 
matrix method for radiative processes is discussed in the OP 
papers by Seaton (1987) and Berringtion et al. (1987). The 
wavefunction W(E) for the (N+l) electron system with total 
spin and orbital angular momenta symmetry S Ln or total angu- 



lar momentum symmetry Jn is expanded in terms of 'frozen' 
jV-electron target ion functions Xi an d vector coupled collision 
electrons 0„ 

<F E (e + ion) = M ^XiQonWi + J] c J / e + ion ^ > W 
> J 

in some specific state S iLim or level y,7r„ index i marking chan- 
nels S iLi(Ji)ni kf£i(S Ln or Jn) with energy k\ of the colliding 
electron described by 9. A is the antisymmetrization operator. 
The second sum expands correlation functions &j as products 
with N+l bound orbital functions that (a) compensate for the 
orthogonality conditions between the continuum and the bound 
orbitals, and (b) represent additional short-range correlation 
that is often of crucial importance in scattering and radiative 
CC calculations for each S Ln. c/s are variational coefficients. 

/^-matrix solutions of coupled equations to total symmetries 
LS states for (N + l)-electron system, 

H N+X V = E¥, (2) 

are discrete bound states W B for E < 0. Seaton (1985) de- 
scribes the matching procedure for inner region and outer re- 
gion wavefunctions at the R-matrix boundary which enables 
calculatons of all required energy levels of the Hamiltonian ma- 
trix. 

For a bound-bound transition from an initial state i to a final 
state j, the line strength, S , is defined as 

s = i < Vj\mWi > t o) 

in atomic units (a.u.), where the dipole operator D is 

D = £V„, (4) 

n 

in the length form, summed over the total number of electrons 
in the ion, and ^ and are the initial and the final wavefunc- 
tions respectively. The oscillator strength fij (/-value) can be 
obtained from line strength S in atomic units ss 

S = (3g i /E ij )f ij . (5) 

where Ejj is the transition energy in rydbergs and g, is the statis- 
tical weight factor of the initial state, with g, = (2S , + 1 )(2L,- + 1 ) 
for a LS multiplet, and (27, + 1) for a fine-structure tran- 
sition. The radiative transition probability, Aj, (Einstein's A- 
coefficient or the radiative decay rate), can be obtained from 
the oscillator strength as 

1 , , Aji(a.u.) 

Ajiia.u.) = -a-E;J n . A^s' 1 ) = '-, (6) 

2 gj T 

where a is the fine-structure constant and t„ = 2.419 lxl0~ 17 s 
is the atomic unit of time. 

The fine-structure components of an LS multiplets can 
be obtained through algebraic transformations of either line 
strength or directly from the oscillator strength (e.g. Nahar 
1995). The line strength, which is independent of the transition 
energy 2?y, is a better choice. As the observed transition ener- 
gies are determined more accurately than the calculated ones, 
use of the former with the line strengths can provide more ac- 
curate fij and Aji. Hence, present fine-structure components are 
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obtained from S values whenever observed energies are avail- 
able, and from /-values when calculated E t j are used. 
The fine-structure line strengths, S jj, are obtained as 

S n = C(J h Jj)S LS /[(25 { + l)(2Lj + + 1)], (7) 

for the allowed transitions (A/ = 0, +1). 5, is the spin which is 
the same as S j. The values of the coefficients C(/„ Jj) can be 
found in Allen (1976). The S jj values satisfy the condition 

Sls=J]Sjj. (8) 
./ 

The fine-structure f-values, ///, can be obtained directly from 
f LS as (Seaton et al 1994): 

fj J (n j S i LjJj,n i SiLiJi) = fisfrjS ,L;,»,S i L i )(2J j + 1)(2L, + 1) 

x W 2 (L j L i J j J i -lS i ), (9) 

where W{LjLiJ jJc 15,) is a Racah coefficient. The above com- 
ponents also satisfy the sum rule 

^(27, + l)fjj(tijS t LjJ j,HiS ,L,J,) = f LS (rijS jLj^jS jLj) x 

(2S, + 1)(2L, + 1). (10) 

The above form for fine structure components is used when 
one or both LS terms of the transition are unobserved or when 
not all the fine structure levels of an LS term are observed, and 
for transitions between high angular momentum states, such 
as higher than H «-» / for which Allen's coefficients are not 
available. 

The lifetime of a level can be computed as 

1 - zk, <"> 

where the sum is the total radiative transition probability for 
level k, i. e. 

giA^l = 2.6774 x 10V (£, - E k ) 3 S E1 (i, k) (12) 

(the observed rate) in the electric dipole case. 

The opacity is a measure of radiation transport through 
matter. The monochromtic opacity is a function of the pho- 
ton frequency and depends mainly on the detailed atomic data: 
bound-bound oscillator strengths and bound-free photoioniza- 
tion cross sections. The quantity of practical interest in stel- 
lar models is the Rosseland mean opacity, K R (T,p), defined as 
(Seaton etal. 1994) 

i r°° i 15 

-= -g(u)du, g ( u )=— u 4 e- u (l-e- u r 2 , (13) 
kr Jo Kv 47T 4 

where g(u) is the Planck weighting function and u = hv/kT. 
Monochromatic opacities k v depends primarily on oscillator 
strengths as 

2 

K v (i -» j) = —Nif(ij)<t> v (14) 
mc 

where N t = ion density in state i, cp v is a profile factor, and on 
photoionization cross sections, o>/, as 



Table 1 . Target terms and relative energies (in Ry) of Fe V in the wave- 
function expansion of Fe IV. 



k v = Nia-pi{v). 



(15) 





Term 




E /Ry 


1 


3d 4 


5 D 


0.000000 


2 


3d 4 


3 H 


0.230232 


3 


3d 4 


3 P 


0.234211 


4 


3d 4 


3 F 


0.244942 


5 


3d 4 


3 G 


0.274991 


6 


3d 4 


>G 


0.333399 


7 


3d 4 


3 D 


0.334726 


8 


3d 4 


i/ 


0.341832 


9 


3d 4 


'5 


0.361166 


10 


3d 4 


l D 


0.421837 


11 


3d 4 


l F 


0.480536 


12 


3d 4 


3 F 


0.567719 


13 


3d 4 


3 P 


0.568463 


14 


3d 4 


l G 


0.649554 


15 


3d 4 


l D 


0.855063 


16 


3d 4 


l S 


1.103820 



We compute the monochromatic and the mean opacities as 
a measure of the differences between the present dataset for the 
bound-bound transitions and the earlier OP data for Fe IV. 



3. Computations 

The R-matrix computations for Fe IV are carried out using a 
wavefunction expansion of 16 terms of target Fe V belonging 
to the ground configuration, 3d 4 , as given in Table 1 ; the en- 
ergies used are the measured values available from the NIST 
compilation at www.nist.gov. The orbital wavefunctions of the 
target were obtained from atomic structure calculations using 
SUPERSTRUCTURE (Eissner et al. 1974), as in the earlier 
works of Bautista and Pradhan (1997) and Nahar et al. (1998). 
The present wavefunction expansion is a subset of the previ- 
ous 3 1 term expansion, but includes singlet target terms which 
were not considered earlier. This enables the computation of 
the large number of doublet (N+l)-electron symmetries. 

In addition to the target expansion in Table 1, it is necessary 
to ensure an accurate representation of short range correlation 
effects that are vital for the bound (N+l)-electron states. The 
second term of the wavefunction expansion, Eq.(l), consists 
of 61 such 'bound channel' correlation functions of Fe IV, as 
described by Bautista and Pradhan (1997). 

The R-matrix calculations span several stages of computa- 
tions (Berrington etal. 1987, 1995). The computations were 
extensive in terms of CPU memory size and time require- 
ment. Hence, they were carried out one or a few symmetries 
of S Ln at a time. The one-electron orbital functions optimised 
through a configuration-interaction type calculation for Fe V 
using SUPERSTRUCTURE are the input to STG1 of the R- 
matrix codes to compute Slater integrals. Radial integrals for 
the partial wave expansion in Eq. 1 are specified for orbitals 
< I < 9 as a basis of NRANG2 = 12 'continuum' functions 
— sufficient for bound electrons with n < 1 1 at a radius of the 
^-matrix box. Computations are carried out for bound symme- 
tries with < L < 9 (both parities) and spin multiplicities, 
2S+1 = 2,4,6. 

The bound energy eigenvalues are obtained from computer 
program known STGB. All bound states of the Hamiltonian 
matrix are scanned for up to n = 1 1 and I = 9 with a fine ef- 
fective quantum number mesh of Av = 0.001. However, due 
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to the computationl complexities stemming from the numerous 
bound states obtained, the code was unable to generate the stan- 
dard R-matrix energy files. A comlementary computer code, 
PLSRAD, was written to sort out the energies directly from the 
eigenvalues of the Hamiltonian matrix. 

Furthermore, identification of the large number of LS bound 
states was a formidable task for this ion. With the excep- 
tion of equivalent-electron states, the energies can be identified 
through quantum defect analysis of Rydberg series: 

z 2 

E = E, — - , v = n - /ui ± i/ 2 (t) (16) 

with series limits E, at the 16 Fe v 'target' states, and the 
weighted percentage contributions of interacting closed chan- 
nels. An LS term is designated according to a possible combi- 
nation of the configuration of the core and the outer electron, 
and assigned an appropriate value of v for the outer electron 
quantum numbers n, I associted with the maximum channel 
percentage contribution. A Rydberg series is identified by sort- 
ing out the series of v for the terms of same configuration but 
with increasing n of the valence electron. Code ELEVID (e.g. 
Nahar 1995) was used for this analysis. Uncertainties in ass- 
ingments of proper configurations are introduced for cases with 
small differences in v for various terms of the same symmetry, 
as the quantum defects for these states are almost the same. 
There is also some uncertainty in the identification of the ad- 
jacent higher or lower terms for such cases. The code STGBB 
is used to compute radiative data for the bound-bound transi- 
tions; the code exploits methods developed by Seaton (1986) to 
evaluate the outer region (> RA) contributions to the radiative 
transition matrix elements. 

The f-, S -, and A- values for the fine structure transitions are 
obtained through algebraic transformations of the LS multiplet 
line strengths using spectroscopiclly observed energies. For the 
cases where the observed energies are not found for all the fine 
structure components of an LS term, the spitting of LS f-values 
is carried out as described above using the code JJTOLS (Nahar 
1995). 

Using the new bound-bound radiative dataset, but retain- 
ing the earlier OP bound-free data (the contribution of bound- 
free is much smaller), we also calculate the monochromatic 
and Rosseland mean and Planck mean opacities as described 
by Seaton etal. (1994). Results are compared with those com- 
puted using the OP data calculated by Sawey and Berrington 
(1992). 



4. Results and Discussion 

In the subections below we describe the results for bound state 
energies of LS terms and fine structure levels levels of Fe IV, 
radiative transitions for all dipole transitions, and the effect of 
the new data on opacities. We note that this entails the calcula- 
tion and identification of the largest number of LS bound terms 
(1,798) computed for any ion under the Opacity or the Iron 
Projects thus far. 



4.1. Energy Levels 

A total of 1,798 bound states of Fe IV of sextet, quartet and 
doublet symmetries were obtained with n <~ 11, I < 9, 
< L < 9 (both even and odd parities), and spin multiplicities 
2S+1 = 2,4,6. All energy levels have been identified and a cor- 
respondence is established with the observed energies where 
available. 

In Table 2, the present bound state energies are compared 
with the measured values (available from the NIST database), 
the OP calculations by Sawey and Berrington (1992), and by 
Bautista and Pradhan (1997). The low-lying energy levels in 
the present dataset are in very good agreement, within 5%, 
with the measured values. However, the higher ones show 
larger differences, the highest one being about 18% for the 
state 3d 4 ( 1 Dl)4p( 2 D°). A reason for such larger uncertainty 
for higher states, especially the odd parity states, could be the 
improper balance for inclusion of configurations which can im- 
prove the lower states considerably, but introduce some 'spec- 
tral repulsion' for the higher levels of the ion. Nonetheless, 
present energies on average are an improvement over the ones 
by Sawey and Berrington (1992) except for the high lying odd 
parity states where their energies show better agreement with 
experiment by a few percent. Present energies are also slightly 
better than those by Bautista and Pradhan (1997). The improve- 
ment in computed energies for the low-lying states is an impor- 
tant fact in the calculation of opacities, as discussed later. 

The 104 observed LS term energies are obtained from sta- 
tistical averaging of the 269 measured fine structure energies 
(NIST). In contrast, the present 1.798 LS states correspond to 
5,331 fine structure levels; with the exception of the observed 
energy levels, the fine structure transition energies are approx- 
imated as explained in the Theory section. 

The complete table of 1,798 bound state energies of Fe IV 
is available electronically. The energies of 104 observed states 
have been replaced by measured values. The bound states have 
been assigned with a prefixed index letter for convenience, as 
shown in sample Table 3 (the letters may not necessarily match 
with those of the NIST table). Following NIST, it is a standard 
practice to designate the lowest 52 bound states with prefixes 
in alphabetical order. For even symmetries, an ascending order 
of alphabet is used, first in the lower case and then in the upper 
case. The opposite order is chosen for the odd parity states, that 
is, descending order of alphabet, first in the lower case and then 
the upper case. States lying above these 52 states do not have 
any such notation; they are recognized by the energies only. 

4.2. LS multiplets 

The LS coupling oscillator strengths are used in various astro- 
physical applications, such as calculating plasma opacities. The 
LS multiplet data with 1,798 bound states of Fe IV comprises 
138,121 dipole allowed oscillator strengths / and coefficients 
of spontaneous decay, the A - values. The full data table is 
available electronically. Since the line strength S does not de- 
pend on the transition energy, the present S values are used to 
recompute the / and A-values (Eq. 5) using observed energies 
whenever available to improve accuracy. 
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Table 2. Comparison of the present energies E c of Fe IV with the ob- 
served binding energies E , and computed energies, E BP (Bautista and 
Pradhan 1997) and E SB (Sawey and Berrington 1992). An * next to a 
state indicates an incomplete multiplet with unobserved fine structure 
components. 

Config SLn E a E c E bp E sb 

(Ry) (Ry) (Ry) (Ry) 



5 

Table 3. Sample of energies in the complete table of energies. 



Config 


SLn 


E(Ry) 


3a5 


aoS 


-4.020U0 


O 7 A / C r\ \ A 7 

3fl4(5L))4fl 


b6S 


-1.39828 


O J A / C r\ \ C J 

3d4(5D)5d 


c6S 


-0.81636 


o J A / C r\ \ f J 

3d4(5D)bd 


d6S 


f\ C A f\f\ 1 

-0.54091 


O J A / C T~\ \ ^ J 

3d4(5D)ld 


e6S 


-0.38568 


3fl4(5D)8fl 


f6S 


-0.28891 


O 7 A / C r\\f\ J 

3d4(5D)9d 


g6S 


f\ r\r\ a AC\ 

-0.22449 


1 J/1 / ~l \ 1 A 7 


h6S 


-0.17945 


o J /i /c r\\ 11 7 

3a4(5D)l Id 


i6S 


-0.14672 


o j /i /c r\\ /i j 

3fl4(5D)4d 


a6P 


-1.55722 


3d4(5D)5d 


inn 

boP 


-0.8/599 


3d4(5L>)6a 


coP 


-0.5/1 15 


3d4{5D)ld 


d6P 


-0.40367 


3d4(5D)8d 


e6P 


-0.30030 


3d4(5D)9d 


f6P 


-0.23216 


3d4(5D)10d 


g6P 


-0.18486 


3d4(5D)lld 


h6P 


-0.15068 



Table 4 presents a sample of the complete electronic output 
file, with the following explanation. The first line provides the 
ion information, [Z-26, number of core electrons = 22]. Note 
that the number of core electrons is one less than that in the 
final ion. This (somewhat confusing) notation is standard in the 
R-matrix calculations and refers to the 'target ion' with which 
another electron forms a bound state. The first two columns, 
[SiLiPi(i) -> SfLfPf(f)], refers to the symmetries of the 
initial and final bound states; the next two columns to the initial 
and final energies, £, and Ef, in Rydbergs, followed by /, gf, 
A-values. The total number of transitions is (138121) is also 
listed at the end of the headings. The A-values are in atomic 
units; they should be divided by 2.419xl0~ 17 to yield A{sec~ l ). 
Designation for a given state in the first two columns consists 
of integers representing the spin multiplicity 'S', which is in 
fact (2S+1); the total angular momentum L is for S, 1 for 
P, etc.; the parity Pi = for an even state and = 1 for and odd 
state; and the fourth integer is the index of the ascending energy 
order computed by STGB for the given symmetry. 

Present results for the Fe IV oscillator strengths are com- 
pared with earlier works in Table 5. Bautista and Pradhan 
(1997) considered transitions among sextets and quartets only 
and hence have relatively smaller set of data compared to the 
OP data. Under the OP, Sawey and Berrington (1992) also pro- 
duced an extensive set of data, a total of 133,650 transitions 
among 1,721 bound states for all three symmetries, sextets, 
quartets, and doublets. Hence, OP data contains about 4,500 
transitions less than the present set. Fawcett (1989) employed 
the semi-empirical Cowan's Hartree-Fock relativistic code with 
optimized Slater parameters to calculate the transition proba- 
bilities of Fe IV. We note that the NIST database of evaluated 
compilation gives radiative decay rates only for a small number 
of forbidden transitions, not allowed ones. Table 5 compares 
present results with these earlier works for some transitions. 
Agreement among all the calculations is good for sextet and 
quartet transitions, within 10%. However, for the doublet tran- 



3d5 6 S 

3d5 4 G 

3d5 A P 

3d5 4 D 

3d5 2 I 

3d5 3 2 D 

3d5 2 2 F 

3d5 4 F 

3d5 2 H 

3d5 2 2 G 

3d5 \ 2 F 

3d5 2 S 

3d5 2 2 D 

3d5 i 2 G 

3d5 2 P 

3d5 l 2 D 
3d4(5D)4s 6 D 



3d4(5D)4s 4 D 

3d4(3P2)4s 4 P 

3d4(3H)4s 4 H 

3d4{3F2)4s 4 F 

3d4(3G)4s 4 G 

3d4(3P2)4s 2 P 

3d4(3H)4s 2 H 

3d4(3F2)4s 2 F 

3d4{3G)4s 2 G 

3d4(3D)4s A D 

3d4{\G2)4s 2 G 

3rf4(l/)4,s 2 I 

3^4(152)4.? 2 S 

3d4(3D)4s 2 D 

3d4(lD2)4s 2 D 

3d4(lF)4s 2 F 

3d4(5D)4p 6 F° 

3d4(5D)4p 6 P° 

3d4(3Pl)4s A P 

3d4(3Fl)4s 4 F 

3d4(5D)4p 4 P° 

3d4(5D)4p 6 D° 

3d4(3Pl)4s 2 P 

3d4(3Fl)4s 2 F 

3d4(5D)4p 4 F° 

3d4(lGl)4s 2 G 
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sitions the present oscillator strengths and those of Sawey and 
Berrington (1992) show similar values, while Fawcett's values 
are somewhat smaller. 

4.3. Fine structure transitions 

Fine structure /- andA-values are presented for 7 12, 120 dipole 
allowed transitions in Fe IV. As in the case of LS multiplets, 
observed energies of fine structure transitions are employed 
whenever available for imporved accuracy. As only a limited 
number of observed levels are available, a smaller set of 2,915 
fine structure transitions corresponding to 620 LS multiplets 
has observed energies. Rest of the transitions, 709,205 in total, 
employ calculated energies; the transition energies are obtained 
as explained in the Theory section. The entire file containing all 
fine structure transitions is available electronically. 

A sample set of /-, S -, and A-values for the dipole allowed 
fine structure transitions in Fe IV is presented in Table 6. The 
first line of each subset corresponds to the LS transition fol- 
lowed by the fine structure components. The letter prefix desig- 
nation of the transitional states in the table correspond to their 
energy positions, as explained above (see Table 3). The ener- 
gies of the initial and final fine structure levels are given in unit 
of crrT 1 , while the transitional energy differences are given in 
Rydbergs. The A-values are in An asterisk (*) below an LS 
state indicates an incomplete set of observed energy levels, and 
an asterisk for the transitional energy indicates that one or both 
the levels are missing from the observed energy set. 

Based on the accuracy of the calculated energies and com- 
parison with other works, LS multiplets should be of better 
accuracy than those in TOPbase and by Bautista and Pradhan 
(1997). The overall accuracy for most of the present f-, S, and 
A-values should be within 20% for strong transitions, with the 
possible exception of weak transitions. Because of low charge 
of the ion, the relativistic effects are expected to be small. 

4.4. Opacities 

We calculate sample radiative opacities of Fe IV for three rea- 
sons. First, as a completeness and consistency check on the 
huge amount of data computed, second, to gauge the potential 
impact on iron opacities, and third, spectral models for deter- 
mination of iron abundances in astrophysical sources. 

Fig. 1 shows the monochromatic opacities in most of the 
wavelength range where Fe IV is a contributor in an astrophys- 
ical plasma in local-thermodynamic-equilibrium (LTE). The 
present results are compared with those using the bound-bound 
OP data by Sawey and Berrington (1992). In order to com- 
pare only the main source of opacity (lines), we have used the 
same bound-free (photoionization) data in both calculations; 
the bound-free opacity is less than 10% of the total. We choose 
the temperature log(T) = 4.5 and the electron density N e = 17.0. 
The Mihalas-Hummer-Dappen equation-of-state (MHD-EOS; 
Mihalas etal. 1998) used in the OP work gives relatively large 
value of the ionization fraction Fe IV/Fe = 0.57; Fe IV dom- 
inates the iron opacity at this temperature and density. Fig. 1 
shows immediately the reason that Fe IV is the dominant ion in 



the UV, between 500-3000 A. It is noteworthy that the opacity 
of Fe IV varies over nearly 10 orders of magnitude, implying 
that it remains a large contributor over a very wide range of 
radiation. 

The two datasets displayed in Fig. 1 differ considerably. 
The main features throughout the oapcity spectrum appear at 
different positions. For example, there is relative displacement 
of nearly 200 Ain the position of the trough which is at ~1500 
Ain the present results, but at ~ 1700 Ausing the OP data. 
Even more striking are the diffrences in the present and the 
OP Rosseland mean opacities, 3.3 and 4.87 cm 2 /g respectively. 
This nearly 50% difference is similar to the one we found ear- 
lier in the Fe II opacities (Nahar and Pradhan 1994) using new 
IP data as opposed to the OP data by Sawey and Berrington 
(1997), which has been shown to be inaccurate for low ioniza- 
tion states of iron Fe I- IV (Nahar and Pradhan 1994, Bautista 
and Pradhan 1997, Cowley and Bautista 2003). We note that 
the calculation of latest stellar opacities under the OP has been 
carried out using the OP data as well as data from other sources 
(M.J. Seaton, private communication). 

Although more extensive investigations are needed to as- 
certain the full range of differences at all temperatures and den- 
sities, we analyze the results somewhat further. The improve- 
ment in computed energies in the present calculations for the 
low-lying states is an important fact in the calculation of opaci- 
ties, since they account for most of the level population; higher 
levels are (a) not populated significantly, and (b) smeared out 
or dissolved due to stark and electron impact broadening as in- 
corporated in MHD-EOS. An examination of the all level pop- 
ulations of Fe IV at T = 3.16 x 10 5 K and N e = 10 17 cm 4 , ob- 
tained from the MHD-EOS shows that there is nearly a sharp 
drop in population of levels as one approaches the highest of 
the 16 ground configuration 3d 5 levels. For example, the pop- 
ulation drops by about 2 orders of magnitude from the ground 
state 3d 5 ( 6 S) at ~ -4 Ryd, to the first of the excited configura- 
tion state 3d 4 4s( 6 D) at ~ -2.8 Ryd. Even higher levels do not 
contribute to the opacity in any significant way. The decrease is 
related to the negative exponential functional form in the mod- 
ified Boltzmann-Saha equations in MHD-EOS. 

5. Conclusions 

Results are presented from a large and extensive calculation 
for Fe IV. All 1,798 states were identified in spectroscopic no- 
tation through quantum defect and channel contribution anal- 
ysis. However, level identifications for highly mixed excited 
states should be treated with caution since often no definitive 
assignment is possible, and in any case not needed for prac- 
tical applications. Present energies agree very well with the 
observed enegies for all low lying terms. However, the uncer- 
tainty increases for excited states with higher energies, espe- 
cially the odd parity states. A small subset of the total num- 
ber of corresponding 138,121 transitions among the LS terms, 
and 712,120 fine structure transitions, was reprocessed wher- 
ever observed energies are available. While the transition prob- 
abilities for the strong transitions should not exceed about 20%, 
the uncertainty in the weak fine structure compoments could be 



S. N. Nahar and A.K. Pradhan: Radiative Data for Fe IV 



7 



much higher, especially if relativistic effects are important and 
deviations from LS coupling are significant. 

The new Fe IV results should be particularly useful for the 
analysis and application to ultraviolet to optical spectra from 
astrophysical and laboratory sources. All data files are elec- 
tronically available from the CDS. 
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Table 4. Oscillator strengths, /, values and radiative decay rates, A for transitions in Fe IV [Z=26 , number of core electrons=22]. The first 
transition, [6 8 1 6 1 1 1], means 6 5 (1) — > 6 P°(l), where the number within parenthesis is the energy position of the state. The energies 
are Rydberg and the A-values are in atomic units, i.e. A(sec~ l = A(a.u.)/2.419 x 1CT 17 . 
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Table 5. Comparison of present /-values for Fe IV with previous works: a - Bautista and Pradhan (1997), b - Sawey and Berrington (1992), c - 
Fawsett(1989). 
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MA) 



Fig. 1. Monochromatic opacities of Fe IV using the present data (top panel) and with the Opacity Project data (lower panel) from 
Sawey and Berrington (1992) at Log T (K) = 4.5 and Log N e (cm -3 ) = 17.0. The major features show a shift in wavelength or 
differences in intensities. 
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Table 6. Sample data for allowed fine structure transitions of Fe IV. 
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189885.00 


5.266E+01 


6 


4 


9.381E-02 


9.759E-01 


3.384E+09 


a 6 D e 


->z 6 P" 


2.8492 


2.2879 


5.613E-01 


30 


18 


1.827E-01 


2.929E+01 


7.704E+08 






128967.000 


190226.00 


1.632E+02 


10 


8 


1.817E-01 


9.762E+00 


5.684E+08 






128541.000 


190226.00 


1.621E+02 


8 


8 


6.554E-02 


2.798E+00 


1.664E+08 






128191.000 


190226.00 


1.612E+02 


6 


8 


1.461E-02 


4.653E-01 


2.813E+07 






128541.000 


190008.00 


1.627E+02 


8 


6 


1.170E-01 


5.011E+00 


3.930E+08 






128191.000 


190008.00 


1.618E+02 


6 


6 


1.120E-01 


3.579E+00 


2.855E+08 






127929.000 


190008.00 


1.611E+02 


4 


6 


5.523E-02 


1.171E+00 


9.464E+07 






128191.000 


189885.00 


1.621E+02 


6 


4 


5.692E-02 


1.822E+00 


2.167E+08 






127929.000 


189885.00 


1.614E+02 


4 


4 


1.286E-01 


2.733E+00 


3.293E+08 






127766.000 


189885.00 


1.610E+02 


2 


4 


1.842E-01 


1.952E+00 


2.371E+08 


a 6 O c 


->z 6 D° 


2.8492 


2.2583 


5.910E-01 


30 


30 


3.085E-01 


4.698E+01 


8.654E+08 






128967.000 


193789.00 


1.543E+02 


10 


10 


2.510E-01 


1.275E+01 


7.035E+08 






128967.000 


193386.00 


1.552E+02 


10 


8 


5.669E-02 


2.897E+00 


1.962E+08 






128541.000 


193789.00 


1.533E+02 


8 


10 


7.178E-02 


2.897E+00 


1.631E+08 






128541.000 


193386.00 


1.542E+02 


8 


8 


1.380E-01 


5.607E+00 


3.872E+08 






128541.000 


192595.00 


1.561E+02 


8 


6 


9.789E-02 


4.025E+00 


3.572E+08 






128541.000 


192595.00 


1.561E+02 


8 


6 


9.789E-02 


4.025E+00 


3.572E+08 






128191.000 


193386.00 


1.534E+02 


6 


8 


1.328E-01 


4.025E+00 


2.825E+08 






128191.000 


192595.00 


1.553E+02 


6 


6 


5.249E-02 


1.610E+00 


1.452E+08 






128191.000 


193271.00 


1.537E+02 


6 


4 


1.238E-01 


3.759E+00 


5.248E+08 






127929.000 


192595.00 


1.546E+02 


4 


6 


1.846E-01 


3.759E+00 


3.432E+08 






127929.000 


193271.00 


1.530E+02 


4 


4 


3.419E-03 


6.891E-02 


9.738E+06 






127929.000 


193120.00 


1.534E+02 


4 


2 


1.210E-01 


2.443E+00 


6.857E+08 






127766.000 


193271.00 


1.527E+02 


2 


4 


2.431E-01 


2.443E+00 


3.478E+08 






127766.000 


193120.00 


1.530E+02 


2 


2 


6.933E-02 


6.985E-01 


1.975E+08 


a 6 D e 


-*z 6 F° 


2.8492 


2.2955 


5.537E-01 


30 


42 


4.142E-01 


6.731E+01 


7.286E+08 






128967.000 


190276.00 


1.631E+02 


10 


12 


3.582E-01 


1.923E+01 


7.483E+08 






128967.000 


189515.00 


1.652E+02 


10 


10 


5.424E-02 


2.949E+00 


1.326E+08 






128967.000 


188904.00 


1.668E+02 


10 


8 


4.319E-03 


2.372E-01 


1.293E+07 






128541.000 


189515.00 


1.640E+02 


8 


10 


3.020E-01 


1.305E+01 


5.992E+08 






128541.000 


188904.00 


1.657E+02 


8 


8 


9.991E-02 


4.359E+00 


2.428E+08 






128541.000 


188428.00 


1.670E+02 


8 


6 


1.246E-02 


5.481E-01 


3.975E+07 






128191.000 


188904.00 


1.647E+02 


6 


8 


2.542E-01 


8.270E+00 


4.687E+08 






128191.000 


188428.00 


1.660E+02 


6 


6 


1.359E-01 


4.456E+00 


3.288E+08 






128191.000 


188086.00 


1.670E+02 


6 


4 


2.333E-02 


7.693E-01 


8.373E+07 






127929.000 


188428.00 


1.653E+02 


4 


6 


2.121E-01 


4.616E+00 


3.451E+08 






127929.000 


188086.00 


1.662E+02 


4 


4 


1.669E-01 


3.654E+00 


4.029E+08 






127929.000 


187878.00 


1.668E+02 


4 


2 


3.240E-02 


7.116E-01 


1.553E+08 






127766.000 


188086.00 


1.658E+02 


2 


4 


1.821E-01 


1.987E+00 


2.209E+08 






127766.000 


187878.00 


1.664E+02 


2 


2 


2.283E-01 


2.500E+00 


5.502E+08 


a A P e 


^z*S" 


3.6983 


1.6735 


2.025E+00 


12 


4 


6.310E-02 


1.122E+00 


6.233E+09 






35253.800 


257503.00 


4.499E+01 


6 


4 


6.311E-02 


5.609E-01 


3.119E+09 






35333.300 


257503.00 


4.501E+01 


4 


4 


6.309E-02 


3.739E-01 


2.077E+09 






35406.600 


257503.00 


4.503E+01 


2 


4 


6.307E-02 


1.870E-01 


1.038E+09 


a Ape _> z 4po 


3.6983 


2.2657 


1.433E+00 


12 


12 


5.704E-02 


1.433E+00 


9.402E+08 






35253.800 


193549.00 


6.317E+01 


6 


6 


4.020E-02 


5.017E-01 


6.720E+08 






35253.800 


191694.00 


6.392E+01 


6 


4 


1.703E-02 


2.150E-01 


4.170E+08 






35333.300 


193549.00 


6.320E+01 


4 


6 


2.583E-02 


2.150E-01 


2.875E+08 






35333.300 


191694.00 


6.395E+01 


4 


4 


7.564E-03 


6.371E-02 


1.234E+08 






35333.300 


191021.00 


6.423E+01 


4 


2 


2.354E-02 


1.991E-01 


7.611E+08 






35406.600 


191694.00 


6.398E+01 


2 


4 


4.726E-02 


1.991E-01 


3.849E+08 






35406.600 


191021.00 


6.426E+01 


2 


2 


9.410E-03 


3.982E-02 


1.520E+08 



